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The motion of simple domain walls and of more complex magnetic textures in the presence of a 
transport current is described by the Landau-Lifshitz-Slonczewski (LLS) equations. Predictions of 
the LLS equations depend sensitively on the ratio between the dimensionless material parameter 
j3 which characterizes non-adiabatic spin-transfer torques and the Gilbert damping parameter a. 
This ratio has been variously estimated to be close to 0, close to 1, and large compared to 1. By 
identifying j3 as the influence of a transport current on a, we derive a concise, explicit and relatively 
simple expression which relates /3 to the band structure and Bloch state lifetimes of a magnetic 
metal. Using this expression we demonstrate that intrinsic spin-orbit interactions lead to intra- 
band contributions to j3 which are often dominant and can be (i) estimated with some confidence 
and (ii) interpreted using the "breathing Fermi surface" model. 

PACS numbers: 



I. INTRODUCTION 

An electric current can influence the magnetic state of 
a ferromagnet by exerting a spin transfer torque(STT) 
on the magnetization^! 3 - This effect occurs whenever 
currents travel through non-collincar magnetic systems 
and is therefore promising for magnetoelectronic appli- 
cations. Indeed, STT's have already been exploited in a 
number of technological devices^ Partly for this reason 
and partly because the quantitative description of order 
parameter manipulation by out-of-equilibrium quasipar- 
ticles poses great theoretical challenges, the study of the 
STT effect has developed into a major research subficld 
of spintronics. 

Spin transfer torques are important in both magnetic 
multilayers, where the magnetization changes abruptly^ 
and in magnetic nanowires, where the magnetization 
changes smoothly £ Theories of the STT in systems with 
smooth magnetic textures identify two different types of 
spin transfer. On one hand, the adiabatic or Slonczewski 3 - 
torque results when quasiparticlc spins follow the under- 
lying magnetic landscape adiabatically. It can be math- 
ematically expressed as (v s ■ V)so, where So stands for 
the magnetization and v s is the "spin velocity" , which 
is proportional to the charge drift velocity, and hence to 
the current and the applied electric field. The micro- 
scopic physics of the Slonczewski spin-torque is thought 
to be well understooc£&£, at leasl£ in systems with weak 
spin-orbit coupling. A simple angular momentum conser- 
vation argument argues that in the absence of spin-orbit 
coupling v s = er s E/eso) where sq is the magnetization, a s 
is the spin conductivity and E is the electric field. How- 
ever, spin-orbit coupling plays an essential role in real 
magnetic materials and hence the validity of this sim- 
ple expression for v s needs to be tested by more rigorous 
calculations. 

The second spin transfer torque in continuous media, 



/3so x (v s ■ V)so, acts in the perpendicular direction and 
is frequently referred to as the non-adiabatic torque^ 
Unfortunately, the name is a misnomer in the present 
context. There are two contributions that have the pro- 
ceeding form. The first is truly non-adiabatic and occurs 
in systems in which the magnetization varies too rapidly 
in space for the spins of the transport electrons to fol- 
low the local magnetization direction as they traverse 
the magnetization texture. For wide domain walls, these 
effects are expected to be smallJ^ The contribution of 
interest in this paper is a dissipative contribution that 
occurs in the adiabatic limit. The adiabatic torque dis- 
cussed above is the reactive contribution in this limit. As 
we discuss below, processes that contribute to magnetic 
damping, whether they derive from spin-orbit coupling 
or spin-dependent scattering, also give a spin-transfer 
torque parameterized by (3 as above. In this paper, we 
follow the common convention and refer to this torque as 
non-adiabatic. However, it should be understood that it 
is a dissipative spin transfer torque that is present in the 
adiabatic limit. 

The non-adiabatic torque plays a key role in current- 
driven domain wall dynamics, where the ratio between 
(3 and the Gilbert parameter a can determine the veloc- 
ity of domain walls under the influence of a transport 
current. There is no consensus on its magnitude of the 
parameter Although there have a few theoretical 

studie a 12 ' 13 ! 14 of the STT in toy models, the relationship 
between toy model STT's and STT's in either transition 
metal fcrromagncts or ferromagnetic semiconductors is 
far from clear. As we will discuss the toy models most 
often studied neglect spin-orbit interactions in the band- 
structure of the perfect crystal, intrinsic spin-orbit inter- 
actions, which can alter STT physics qualitatively. 

The main objectives of this paper are (i) to shed new 
light on the physical meaning of the non-adiabatic STT 
by relating it to the change in magnetization damping 
due to a transport current, (ii) to derive a concise for- 
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mula that can be used to evaluate f3 in real materials 
from first principles and (iii) to demonstrate that a and 
(3 have the same qualitative dependence on disorder (or 
temperature), even though their ratio depends on the 
details of the band structure. As a byproduct of our the- 
oretical study, we find that the expression for v s in terms 
of the spin conductivity may not always be accurate in 
materials with strong spin-orbit coupling. 

We begin in Section II by reviewing and expanding 
on microscopic theories of a, (3 and v s . In short, our 
microscopic approach quantifies how the micromagnctic 
energy of an inhomogeneous ferromagnet is altered in 
response to external rf fields and dc transport currents 
which drive the magnetization direction away from lo- 
cal equilibrium. These effects are captured by the spin 
transfer torques, damping torques, and effective magnetic 
fields that appear in the LLS equation. By relating mag- 
netization dynamics to effective magnetic fields, we de- 
rive explicit expressions for a, (3 and v s in terms of mi- 
croscopic parameters. Important contributions to these 
materials parameters can be understood in clear physical 
terms using the breathing Fermi surface modelJ^ Read- 
ers mainly interested in a qualitative explanation for our 
findings may skip directly to Section VIII where we dis- 
cuss of our main results in that framework. Regardless 
of the approach, the non-adiabatic STT can be under- 
stood as the change in the Gilbert damping contribution 
to magnetization dynamics when the Fermi sea quasi- 
particlc distribution function is altered by the transport 
electric field. The outcome of this insight is a concise an- 
alytical formula for f3 which is simple enough that it can 
be conveniently combined with first-principles electronic 
structure calculations to predict /3-values in particular 
materials^ 

In Sections III, IV and V we apply our expression for (3 
to model ferromagnets. In Section III we perform a nec- 
essary reality check by applying our theory of (3 to the 
parabolic band Stoner ferromagnet, the only model for 
which more rigorous fully microscopic calculations^^ 
of (3 have been completed. Section IV is devoted to 
the study of a two-dimensional electron-gas ferromag- 
net with Rashba spin-orbit interactions. Studies of this 
model provide a qualitative indication of the influence 
of intrinsic spin-orbit interactions on the non-adiabatic 
STT. We find that, as in the microscopic theor y 17 i 18 
for a, spin-orbit interactions induce intra-band contri- 
butions to P which are proportional to the quasiparticlc 
lifetimes. These considerations carry over to the more 
sophisticated 4-band spherical model that we analyze in 
Section V; there our calculation is tailored to (Ga,Mn)As. 
We show that intra-band (conductivity-like) contribu- 
tions arc prominent in the 4-band model for experimen- 
tally relevant scattering rates. 

Section VI discusses the phenomenologically important 
a I ' (3 ratio for real materials. Using our analytical results 
derived in Section II (or Section VIII) we are able to re- 
produce and extrapolate trends expected from toy mod- 
els which indicate that a/ (3 should vary across materials 



in approximately the same way as the ratio between the 
itinerant spin density and the total spin density. We also 
suggest that a and (3 may have the opposite signs in sys- 
tems with both hole-like and electron-like carriers. We 
present concrete results for (Ga,Mn)As, where we obtain 
a I ' (3 ~ 0.1. This is reasonable in view of the weak spin 
polarization and the strong spin-orbit coupling of valence 
band holes in this material. 

Section VII describes the generalization of the torque- 
correlation formula employed in ab-initio calculations of 
the Gilbert damping to the case of the non-adiabatic 
spin-transfer torque. The torque correlation formula in- 
corporates scattering of quasiparticles simply by intro- 
ducing a phcnomcnological lifetime for the Bloch states 
and assumes that the most important electronic transi- 
tions occur between states near the Fermi surface in the 
same band. Our ability to make quantitative predictions 
based on this formula is limited mainly by an incomplete 
understanding of Bloch state scattering processes in real 
ferromagnetic materials. These simplifications give rise 
to ambiguities and inaccuracies that we dissect in Section 
VII. Our assessment indicates that the torque correlation 
formula for (3 is most accurate at low disorder and rela- 
tively weak spin-orbit interactions. 

Section VIII restates and complements the effective 
field calculation explained in Section II. Within the adi- 
abatic approximation, the instantaneous energy of a fer- 
romagnet may be written in terms of the instantaneous 
occupation factors of quasiparticle states. We determine 
the effect of the external perturbations on the occupation 
factors by combining the relaxation time approximation 
and the master equation. In this way we recover the re- 
sults of Section II and are able to interpret the intra-band 
contributions to (3 in terms of a generalized breathing 
Fermi surface picture. 

Section IX contains a brief summary which concludes 
this work. 



II. MICROSCOPIC THEORY OF a, (3 AND v s 

The Gilbert damping parameter a, the non-adiabatic 
spin transfer torque coefficient (3 and the "spin velocity" 
v s appear in the generalized Landau-Lifshitz-Gilbert ex- 
pression for collective magnetization dynamics of a fer- 
romagnet under the influence of an electric current: 

(d t + v s • V) Cl - Cl x H eS = x d t ti -j3Clx (v s • V)fi. 

(1) 

In Eq. ([T]) Ties is an effective magnetic field which we 
elaborate on below and Cl = So/so — (Ct x , Cl y , 1 — + 
D,y)/2) is the direction of the magnetization^ Eq. ([T]) de- 
scribes the slow dynamics of smooth magnetization tex- 
tures in the presence of a weak electric field which induces 
transport currents. It explicitly neglects the dynamics of 
the magnetization magnitude which is implicitly assumed 
to be negligible. For small deviations from the easy di- 
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rection (which we take to be the z-direction), it reads 
H e s,x = (d t + v s • V) n y + (ad t + /3v 8 • V) Q x 

n eSlV = -(d t + v s - v) n x + (ad t + /3v s • v) n y (2) 

The gyromagnetic ratio has been absorbed into the units 
of the field 7i e ff so that this quantity has inverse time 
units. We set ft = 1 throughout. 

In this section we relate the a, (5 and v s parameters 
to microscopic features of the ferromagnet by consider- 
ing the transverse total spin response function. For a 
technically more accessible (yet less rigorous) theory of 



a and p we refer to Section VIII. The transverse spin re- 
sponse function studied here describes the change in the 
micromagnetic energy due to the departure of the magne- 
tization away from its equilibrium direction, where equi- 
librium is characterized by the absence of currents and 
external rf fields. This change in energy defines an ef- 
fective magnetic field which may then be identified with 
Eq. thereby allowing us to microscopically determine 
a,P and v s . To first order in frequency uj, wave vector q 
and electric field, the transverse spin response function 
is given by 



So&a = ^ Xa.bH. 



cxt.b 



Ti-cxt.b 



(3) 



where a,b € {x,y}, 7i GX t is the external magnetic field with frequency to and wave vector q, Sq = sqV is the total 
spin of the ferromagnet (V is the sample volume), and x is the transverse spin-spin response function in the presence 
of a uniform time-independent electric field: 

Xa, 6 (q,a;;v 8 ) = i^ dt J dr exp(ia;t - iq ■ r)([S a (r, t), S b (0, 0)]). (4) 
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In Eq. (3|), x (0) = X(q = 0,w = 0;E = O) de- 
scribes the spin response to a constant, uniform ex- 
ternal magnetic field in absence of a current, x = 
lim^^o x( c t — 0; u \ E = 0)/w characterizes the spin re- 
sponse to a timc-dcpcndcnt, uniform external mag- 
netic field in absence of a current, and x^ = 
lim g) „ s ^o x(q, oj = 0;E)/q • v s represents the spin re- 
sponse to a constant, non- uniform external magnetic field 
combined with a constant, uniform electric field E . Note 
that first order terms in q are allowed by symmetry in 
presence of an electric field. In addition, () is a ther- 
mal and quantum mechanical average over states that 
describe a uniformly magnetized, current carrying ferro- 
magnet. 

The approach underlying Eq. ([3]) comprises a linear 
response theory with respect to an inhomogeneous mag- 
netic field followed by a linear response theory with re- 
spect to an electric field. Alternatively, one may treat the 
electric and magnetic perturbations on an equal footing 
without predetermined ordering; for further considera- 
tions on this matter we refer to Appendix A. 

In the following we emulate and appropriately gen- 
eralize a procedure outlined elsewhereJi First, we rec- 
ognize that in the static limit and in absence of a cur- 
rent the transverse magnetization responds to the exter- 
nal magnetic field by adjusting its orientation to min- 
imize the total energy including the internal energy 
Eint and the energy due to coupling with the exter- 
nal magnetic field, E ext = —So^l ■ TL cx t- It follows 
that xTb = Sl^E^/dtiadClb)- 1 and thus H int , a = 



-(l/S )dE int /dCl a = -Solx^]-]^, where H int is the 
internal energy contribution to the effective magnetic 
field. Multiplying Eq. © on the left by [x (0) ] _1 and 
using H c ff = Hint + Ti-ext we obtain a formal equation for 
'Hen'- 



where 



We ff ,a = £[4> + 45(V S -V) 



(5) 



(6) 



Identifying of Eqs. ([5]) and ^ results in concise micro- 
scopic expressions for a and (3 and v s : 
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In the third line of Eq. ([7]) we have combined the second 
line of Eq. (J6]) with x^ = xl ( v s ■ q)- 

When applying Eq. ((?]) to realistic conducting fer- 
romagnets, one must invariably adopt a self-consistent 
mean-field (Stoner) theory description of the magnetic 
state derived within a spin-density-functional theory 
(SDFT) framework32i In SDFT the transverse spin 
response function is expressed in terms of Kohn-Sham 
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quasiparticlc response to both external and induced mag- 
netic fields; this allows us to transform^ Eq. (|7|) into 

1 Im[x? P _(q = 0,^,E = 0)] 

a = — hm 

bo oj 

1 Im[x? P _(q,^ = 0,E)] 
p = -— hm 

o v s ,q->0 q • v s 

vs-q = -^Rc[x? P _(q,^ = 0,E)], (8) 

Do 

where we have used^ xfl = $a,bSo/A and 

J z — ' £j — e,j — UJ — 177 

(i|5+Ao(r)e^ r |i)(i|5-A (r)e-^ r |i) 

(9) 

is the quasiparticle response to changes in the direction 
of the exchange-correlation effective magnetic field.— To 
estimate (3 this response function should be evaluated 
in the presence of an electric current. In the derivation 
of Eq. IS} we have made use of the fact that Xx,x and 

Xx% are purely imaginary, whereas Xx}y is purely real; 
this can be verified mathematically through S ± = S x ± 
iS y . Physically, "Im" and "Re" indicate that the Gilbert 
damping and the non-adiabatic STT are dissipative while 
the adiabatic STT is reactive. Furthermore, in the third 
line it is implicit that we expand Rc[x^ p ] to first order 
in q and E. 

In Eq. ©, S ± is the spin-rising/lowering operator, 
e, and /j are the Kohn-Sham cigcnstates, eigenenergies 
and Fermi factors in presence of spin-dependent disorder, 
and Ao(r) is the difference in the magnetic ground state 
between the majority spin and minority spin exchange- 
correlation potential - the spin-splitting potential. This 
quantity is always spatially inhomogcncous at the atomic 
scale and is typically larger in atomic regions than in 
interstitial regions. Although the spatial dependence of 
Ao(r) plays a crucial role in realistic ferromagncts, we 
replace it by a phenomenological constant Aq in the toy 
models we discuss below. 

Our expression of v s in terms of the transverse spin 
response function may be unfamiliar to readers familiar 
with the argument given in the introduction of this pa- 
per in which v s is determined by the divergence in spin 



current. This argument is based on the assumption that 
the (transverse) angular momentum lost by spin polar- 
ized electrons traversing an inhomogcncous fcrromagnct 
is transferred to the magnetization. However, this as- 
sumption fails when spin angular momentum is not con- 
served as it is not in the presence of spin-orbit coupling. 
In general, part of the transverse spin polarization lost 
by the current carrying quasiparticles is transferred to 
the lattice rather than to collective magnetic degrees of 
freedom^ when spin-orbit interactions are present. It is 
often stated that the physics of spin non-conservation is 
captured by the non-adiabatic STT; however, the non- 
adiabatic STT per se is limited to dissipative processes 
and cannot describe the changes in the reactive spin 
torque due to spin- flip events. Our expression in terms 
of the transverse spin response function does not rely on 
spin conservation, and while it agrees with the conven- 
tional picture" in simplest cases (see below), it departs 
from it when e.g. intrinsic spin-orbit interactions are 
strong. 

In this paper we incorporate the influence of an electric 
field by simply shifting the Kohn-Sham orbital occupa- 
tion factors to account for the energy deviation of the 
distribution function in a drifting Fermi sea: 

fi ~ /(°>(e, + V t ) ~ /(°)(e 4 ) + VidfW/dei (10) 



where Vi is the effective energy shift for the i-ih eigenen- 
crgy due to acceleration between scattering events by an 
electric field and /(°) is the equilibrium Fermi factor. 
This approximation to the steady-state induced by an 
external electric field is known to be reasonably accurate 
in many circumstances, for example in theories of electri- 
cal transport properties, and it can be used^ to provide 
a microscopic derivation of the adiabatic spin-transfer 
torque. As we discuss below, this ansatz provides a result 
for (3 which is sufficiently simple that it can be combined 
with realistic ab initio electronic structure calculations 
to estimate (3 values in particular magnetic metals. We 
support this ansatz by demonstrating that it agrees with 
full non-linear response calculations in the case of toy 
models for which results are available. 

Using the Cauchy identity, l/(x — if]) = 1/x + iit8(x), 
and df(°' /de ~ — d(e) we obtain 



Im[x? P _] * ^[ W -^ i ]|(i|5 + Ao(r)e^ r K)|^(e i -e F )^-e F ) 



Rc[x? P _] * -jElO'l 5 ' +A °(rK q - r 



2 VjSfa - ep) - Vi8(e t - e F ) 



(11) 



where we have defined the difference in transport devia- tion energies by 



(12) 
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In the first line of Eq. (fTTj) . the two terms within the 
square brackets correspond to the energy of particle- 
hole excitations induced by radio frequency magnetic 
and static electric fields, respectively. The imaginary 

I 



a = 
P = 
v s • q = 



where we have assumed a uniform precession mode for 
the Gilbert damping. 

Eq. KT3J) and Eq. (TT]) identify the non-adiabatic STT 
as a correction to the Gilbert damping in the presence 
of an electric current; in other words, the magnetiza- 
tion damping at finite current is given by the sum of 
the Gilbert damping and the non-adiabatic STT. We feel 
that this simple interpretation of the non-adiabatic spin- 
transfer torque has not received sufficient emphasis in the 
literature. 

Strictly speaking the influence of a transport current 
on magnetization dynamics should be calculated by con- 
sidering non-linear response of transverse spin to both 
effective magnetic fields and the external electric field 
which drives the transport current. Our approach, in 
which we simply alter the occupation probabilities which 
appear in the transverse spin response function is admit- 
tedly somewhat heuristic. We demonstrate below that 
it gives approximately the same result as the complete 
calculation for the case of the very simplistic model for 
which that complete calculation has been carried out. 

In Eq. (fT3")) . the cigenstates indexed by i are not Bloch 
states of a periodic potential but instead the cigenstates 
of the Hamiltonian that includes all of the static dis- 
order. Although Eq. Q13p provides compact expressions 
valid for arbitrary metallic ferromagnets, its practical- 
ity is hampered by the fact that the characterization of 
disorder is normally not precise enough to permit a reli- 

I 



part selects scattering processes that relax the spin of 
the particle-hole pairs mediated either by phonons or by 
magnetic impurities.— Substituting Eq. (fTTj) into Eq. ([8|) 
we can readily extract a, (3 and v s : 



- e F )S(ej - e F ) 



able solution of the Kohn-Sham equations with arbitrary 
impurities. An approximate yet more tractable treat- 
ment of disorder consists of the following steps: (i) re- 
place the actual eigenstates of the disordered system by 
Bloch eigenstates corresponding to a pure crystal, e.g. 
\i) — > |k, a), where k is the crystal momentum and a 
is the band index of the perfect crystal; (ii) switch Vi to 
V a = Tk, a Vk, a ■ eE. where r is the Bloch state lifetime and 
Vk,a = <9ek,a/<9k is the quasiparticlc group velocity, (iii) 
substitute the <5(ek, a — £f) spectral function of a Bloch 
state by a broadened spectral function evaluated at the 
Fermi energy: <5(e k , a - £f) — ► A a (e F ,k)/(2n), where 

A a (e F , k) = ^ -g- (14) 

(cf - e Ka f + -S* 

and r a ,k = 1/ T a,k is the inverse of the quasiparticle 
lifetime. This minimal prescription can be augmented 
by introducing impurity vertex corrections in one of the 
spin-flip operators, which restores an exact treatment of 
disorder in the limit of dilute impurities. This task is 
for the most part beyond the scope of this paper (see 
next section, however). The expression for a in Eq. (|13p 
has already been discussed in a previous paper^i hence 
from here on we shall concentrate on the expression for 
f3 which now reads 



— |015+A (r)|i)| 2 5(e, - e F )S(e 3 e F ) 

i,j 

-1- £ \ti\S+A (ry**\i)f V ^ ' £F) ~ ViS{Ci 
zoo 'rr £j — €j 
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/3 (0) = lim^E / I (k + q, b\S+ A (r)|k, a) \*A a (e F , k)A b (e F ,k + q) ( Vk +^ T W V ^ T ^ \ eE (15) 

where we have used J2 k — » V J d D k/(2Tr) D = V J k with D as the dimensionality V as the volume and 

q . Vs = ±y f I (k + q, b\S+ A (r)|k, a) | 2 e ^+^+^ S ^ ^] ev k , a T k , Q <5(6 F - e k , a ) ^ 

2S0 ^ Jk £k,a — ek+q,6 
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In Eq. (fT5|) the superscript "0" is to remind of the absence of impurity vertex corrections; . In addition, we recall that 
So = So/ V is the magnetization of the ferromagnet and \ak) is a band eigenstate of the ferromagnet without disorder. 
It is straightforward to show that Eq. (|16j) reduces to the usual expression v s = a s E/(eso) for vanishing intrinsic 
spin-orbit coupling. However, we find that in presence of spin-orbit interaction Eq. p6[) is no longer connected to 
the spin conductivity. Determining the precise way in which Eq. (|16[) departs from the conventional formula in real 
materials is an open problem that may have fundamental and practical repercussions. Expanding the integrand in 
Eq. (115|) to first order in q and rearranging the result we arrive at 



?(o) _ 



8irs 



Airs 



V / [|(a,k|5 + A (r)|fe,k)| 2 + |(a,k|S-A (r)|&,k)| 2 ] A a (e F , k)A' b (e F> k)(v k , a • eE)(v k . b • q)r a 

~ E / Re [(fc,k|5-A (r)|a,k)(a,k|5+A (r)q.9 k |6,k) 



(S+ «-> S-)]A a (e F ,k)A b (e F ,k)(v Ka • eE)r a 



(17) 



(a) 
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ent simplified models of ferromagnets. For a simpler-to- 
implement approximate version of Eq. (|15|) or Eq. (|17p 
we refer to Section VI. 

III. NON-ADIABATIC STT FOR THE 
PARABOLIC TWO-BAND FERROMAGNET 



(b) 



eV, 



b, k+q; co+eV ab 




a,k;con 



The model described in this section bears little resem- 
blance to any real ferromagnet. Yet, it is the only model 
in which rigorous microscopic results for (3 are presently 
available, thus providing a valuable test bed for Eq. (|17|) . 
The mean-field Hamiltonian for itinerant carriers in a 
two-band Stoner model with parabolic bands is simply 



FIG. 1: Feynman diagrams for (a) a and (b) /3(q • v s ), the 
latter with a heuristic consideration of the electric field (for a 
more rigorous treatment see Appendix A). Solid lines corre- 
spond to Green's functions of the band quasiparticles in the 
Born approximation, dashed lines stand for the magnon of fre- 
quency uj and wavevector q, uj n is the Matsubara frequency 
and eV a ,b is the difference in the transport deviation energies. 



where A'(e F ,k) = 2(e F - e k , a )r a / [(e F - e k , a ) 2 + T 2 Jl] 2 
stands for the derivative of the spectral function and we 
have neglected <9r/<9k. Eq. (T7]) (or Eq. (|T5])) is the cen- 
tral result of this work and it provides a gateway to eval- 
uate the non-adiabatic STT in materials with complex 
band structures^ for a diagrammatic interpretation see 
Fig. ([1]) . An alternative formula with a similar aspiration 
has been proposed recently^ yet that formula ignores 
intrinsic spin-orbit interactions and relies on a detailed 
knowledge of the disorder scattering mechanisms. In the 
following three sections we apply Eq. (fT7|) to three differ- 



H (k) = ^_ AqS z 
2m 



(18) 



where A is the exchange field and S* b — (5 a h sgn(a). 
In this model the eigenstates have no momentum depen- 
dence and hence Eq. (fT7)) simplifies to 

(v s .q)/3(°) = / A a (e F ,k)A'_ a (e F ,k) 

2tts ^ J k 



k q k • eE 



m m 



-Tk,a, 



(19) 



where a = +(— ) for majority (minority) spins, v k = |- = 
k/m, and S ± = S x ± iS y with S^ b = S Ulb . Also, from 
here on repeated indexes will imply a sum. Taking Ao < 
E F and Ao >> 1/r, the momentum integral in Eq. p^]) 
is performed in the complex energy plane using a keyhole 
contour around the branch cut that stems from the 3D 
density of states: 
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(v B .q)/J<°> = - 9 ° T' q / ^)A a (£ F , a - e )A'_ Q ( eF ._ a - e )er k . a 



eEq 

sgn(a)j/ a e Fja T a r_ a 



6mA .so 
eE ■ q 



("-T r T7l ~ ^in7t) (2°) 



where £F, a = £f + sgn(a)Ao, f a is the spin-dependent 
density of states at the Fermi surface, n a = 2v a e Fta /3 
is the corresponding number density, and 7 a = T a /2. 
The factor 1/3 on the first line of Eq. (j2"0|) comes from 
the angular integration. In the second line of Eq. (|2"0|) 
we have neglected a term that is smaller than the one 
retained by a factor of Aq/(12c|); such extra term (which 
would have been absent in a two-dimensional version of 
the model) appears to be missing in previous work. 13 ' 

The simplicity of this model enables a partial incorpo- 
ration of impurity vertex corrections. By adding to /3^°- ) 



the contribution from the leading order vertex correction 
(P^ 1 '), we shall recover the results obtained previously 
for this model by a full calculation of the transverse spin 
response function. As it turns out, /?W is qualitatively 
important because it ensures that only spin-dependent 
impurities contribute to the non-adiabatic STT in the ab- 
sence of an intrinsic spin-orbit interaction. In Appendix 
B we derive the following result: 



(v 8 • q)/3 



(i) 



Airsn 



u l Rc 



Jk,k' 



D aM D b,b' D b'a' D a',a 



A a (e F ,k) 

(CF - £fe',a') 



A b (e F ,k + q) 



(ep — £k'+q,6') 



V b , a 



A b ,{e Fl k' + q) 

(CF — e k+q,&) 



V v , a 



(21) 



where u l = riiW~ 



0, x, y 7 z), Hi is the density of scatterers, u>, is the Fourier transform of the scattering potential 



and the overline denotes an average over different disorder configurations^ Also, V a ,b 
Expanding Eq. (f2"Tj) to first order in q, we arrive at 



(v, ■ q)/3< 



^a(e F ,k) 



k,k' £ F - £k',a 



A'_ a (e F ,k) A_ Q (e F ,k') 



e F - e k ' 



(CF - £k,-a) 



= (T6«k+q,6 - T- a l> k ,a) ' eE. 

k • q k • eE 



"Tk,a 



(22) 



In the derivation of Eq. (|22j) we have used S 



S x ± iS' !; and assumed that u x 



u x,y , so that 



u l Rc 



QX Qi QX Qi 

a,b D b,b lD b',a'°a'a 



u z ) S a . a '5b.b'5a,-b- In addition, we have used J* k k , -F(|k|, |k = 0. The first 



term inside the square brackets of Eq. (|22[) can be ignored in the weak disorder regime because its contribution is 
linear in the scattering rate, as opposed to the second term, which contributes at zeroth order. Then, 



(v. ■ q)/? (1) 



A 2 

TTSq 



A °V-^) 2eE ' q 



k,k' 



^ a (6F,k)A_ Q (6 F ,k / ) kq keE 

- £k',a)(£F - £k,-a) 2 TTl 



-Tk.o 



7TS 



3m 



j j I ( \ ( /N^a(fF,a - e)A_ a (e F ,-a - e') 

aeae u(e)u(e ) ; — — er„ 



( £ f - e' a )(e F - e-a) g 



/ o ^ eE-q 

~ — 7r(u — u )- — t sign(a)n a r a i/_ a 

2?nAoSo 

Combining this with Eq. (f2"0)) . we get 

(v s • q)/3 ~ (v s ■ q)/3(°' + (v„ • q)/3 (1) 



(23) 



eE • q 

2msoAo 
eEq 

7T — 

ms A 



[n T r T (w 2 ^ + u x ' v i"[) ~ n^i (u z vy + u x ' v v{)} 



(24) 



where we have used j a — n Uu° + u z )v a + 2u x,v v-^ . In ically, whereupon Eq. (|24p agrees with the results pub- 



this model it is simple to solve Eq. ([16]) for v s analyt 
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lished by other authors in Rcfs.[ [l3lfl4| from full non- 
linear response function calculations. However, we reit- 
erate that in order to reach such agreement we had to 
neglect a term of order Ajj/ep in Eq. ([2TJ]) . This extra 
term is insignificant in all but nearly half metallic ferro- 
magnets. 

IV. NON-ADIABATIC STT FOR A 
MAGNETIZED TWO-DIMENSIONAL 
ELECTRON GAS 

The model studied in the previous section misses the 
intrinsic spin-orbit interaction that is inevitably present 
in the band structure of actual ferromagncts. Further- 
more, since intrinsic spin-orbit interaction is instrumental 
for the Gilbert damping at low temperatures, a similarly 
prominent role may be expected in regards to the non- 
adiabatic spin transfer torque. Hence, the present section 
is devoted to investigate the relatively unexplore d 26 ' 27 ef- 
fect of intrinsic spin-orbit interaction on (3. The minimal 



model for this enterprise is the two-dimensional electron- 
gas ferromagnet with Rashba spin-orbit interaction, rep- 
resented by 

= — -b-S, (25) 
2m y ' 

where b = (Xk y , — Xk x , Ao), A is the Rashba spin-orbit 
coupling strength and Ao is the exchange field. 

The eigenspinors of this model are | + ,k) = 
(cos(0/2),-£exp(i0)sin(0/2)) and |-,k) = 
(sin(0/2), i exp(i<j)) cos(0/2)), where the spinor an- 
gles are defined through cos 9 = Ao/i/A 2 fc 2 + A 2 , 
and t&ncf) = k y /k x . The corresponding eigenen- 

ergies are E k ± = k 2 /(2m) =p V + ^ 2 ^ 2 - 
Therefore, the band velocities are given by 

v k± = k(l/m=F A 2 / V /A 2 I 2 TA 2 ) = k/m±. Dis- 
regarding the vertex corrections, the non-adiabatic 
spin-torque of this model may be evaluated analytically 
starting from Eq. (fTT|) . We find that (sec Appendix C): 



(v s • q)/3 



(0) 



A 2 eE- 
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TO 2 


8irsa 




4m + m_ 


A 2 eE- 
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' 1 m 2 \ 2 k 


Sttsq 


2m\ b 2 



b 2 J b 2 

5m + 



1 \ 2 k 2 A 2 



X 2 
b 2 



1 m 2 \ 2 kl 
2mi b 2 



(5m_ A 2 , 
m b 2 



(26) 



where b = ^/A 2 fcp + AQ (fcp = ^/2mep), and 5m± = 
in — m± . As we explain in the Appendix, Eq. (|26[) ap- 
plies for A/cf, Ao, 1/t << £f; for a more general analysis, 
Eq. (fT7|) must be solved numerically (e.g. see Fig. ((2])). 
Eq. (|26|) reveals that intrinsic spin-orbit interaction en- 
ables intra-band contributions to /3, whose signature is 
the 0(t 2 ) dependence on the second line. In contrast, 
the inter-band contributions appear as 0(t ). Since v s 
itself is linear in the scattering time, it follows that (3 
is proportional to the electrical conductivity in the clean 
regime and the resistivity in the disordered regime, much 
like the Gilbert damping a. We expect this qualitative 
feature to be model-independent and applicable to real 
ferromagnets. 



V. NON-ADIABATIC STT FOR (Ga,Mn)As 



In this section we shall apply Eq. (jT7|) to a more sophis- 
ticated model which provides a reasonable description of 
(III,Mn)V magnetic semiconductors^ Since the orbitals 
at the Fermi energy are very similar to the states near 
the top of the valence band of the host (III,V) semicon- 
ductor, the electronic structure of (III,Mn)V ferromag- 
nets is remarkably simple. Using a p-d mean field theory 
model for the ferromagnetic ground state and a four-band 



spherical model for the host semiconductor band struc- 
ture, Gai-^Mn^As may be described by 



H (k) 



1 

2m 



71 



:72 



2 73 (k-S) 



A S 2 



(27) 

where S is the spin operator projected onto the 3—3/2 to- 
tal angular momentum subspace at the top of the valence 
band and {71 = 6.98,72 = 73 = 2.5} are the Luttingcr 
parameters for the spherical approximation to the valence 
bands of GaAs. In addition, Ao = JpdsiVMn = JpdSo is 
the exchange field, J pc j = 55 meVnm 3 is the p-d exchange 
coupling, s = 5/2 is the spin of Mn ions, Nmu = 4s/a 3 
is the density of Mn ions and a = 0.565 nm is the lattice 
constant of GaAs. We solve Eq. (f2"T)l numerically and 
input the outcome in Eqs. (fTB|) . fll7[l . 

The results are summarized in Fig. We find that 
the intra-band contribution dominates as a consequence 
of the strong intrinsic spin-orbit interaction, much like for 
the Gilbert damping^. Incidentally, (3 barely changes 
regardless of whether the applied electric field is along 
the easy axis of the magnetization or perpendicular to it. 
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FIG. 2: M2DEG: inter-band contribution, intra-band con- 
tribution and the total non-adiabatic STT for a magnetized 
two-dimensional electron gas (M2DEG). In this figure the ex- 
change field dominates over the spin-orbit splitting. At higher 
disorder the inter-band part (proportional to resistivity) dom- 
inates, while at low disorder the inter-band part (proportional 
to conductivity) overtakes. For simplicity, the scattering time 
r is taken to be the same for all sub-bands. 
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FIG. 3: GaMnAs: /3^ ' for E perpendicular to the easy axis 
of magnetization (£). x and p are the Mn fraction and the 
hole density, respectively. The intra-band contribution is con- 
siderably larger than the inter-band contribution, due to the 
strong intrinsic spin-orbit interaction. Since the 4-band model 
typically overestimates the influence of intrinsic spin-orbit in- 
teraction, it is likely that the dominion of intra-band con- 
tributions be reduced in the more accurate 6-band model. 
By evaluating j3 for E||z (not shown) we infer that it does 
not depend significantly on the relative direction between the 
magnetic easy axis and the electric field. 



VI. a//3 IN REAL MATERIALS 

The preceding three sections have been focused on test- 
ing and analyzing Eq. (fTT)) for specific models of ferro- 



magnets. In this section wc return to more general con- 
siderations and survey the phenomenologically important 
quantitative relationship between a and f3 in realistic fer- 
romagnets, which always have intrinsic spin-orbit inter- 
actions. We begin by recollecting the expression for the 
Gilbert damping coefficient derived elsewhere^ 



87TS ^—f 
a.b 



f k\S+ Ao\a, k) | 2 A a (e F , k)A b (e Fl k) 

Jk 



(28) 

where we have ignored disorder vertex corrections. This 
expression is to be compared with Eq. (|15[) : for peda- 
gogical purposes we discuss intra-band and inter-band 
contributions separately. 

Starting from Eq. (|15[) and expanding the integrand to 
first order in q wc obtain 



8ttsq 



|(a,k|5+A |a,k)| 2 A a (e F ,k) 2 

eT a q l d ki vl a E j 



q v s 



(29) 



where we have neglected the momentum dependence of 
the scattering lifetime and a sum over repeated indices 
is implied. Remarkably, only matrix elements that arc 
diagonal in momentum space contribute to Antra ; the 
implications of this will be highlighted in the next section. 
Recognizing that dk j v l ka = (l/m)^, where (l/m) a is 
the inverse effective mass tensor corresponding to band 
a, Eq. (|29|) can be rewritten as 



Antra — 



where 



8tts 



Jk 



Ka,k|5 + A |a,k)| 2 A a (e F ,k) 



2 q • Vd,a 



q v s 



v* = er.tm- 1 )^ 



(30) 



(31) 



is the "drift velocity" corresponding to the quasiparticles 
in band o. For Galilean invariant systems^ «d,a = v s 
for any (k, a) and consequently Antra = centra- At first 
glance, it might appear that v s , which (at least in ab- 
sence of spin-orbit interaction) is determined by the spin 
current, must be different than Vd. a - However, recall that 
v s is determined by the ratio of the spin current to the 
magnetization. If the same electrons contribute to the 
transport as to the magnetization, v s = «d, a provided 
the scattering rates and the masses are the same for all 
states. These conditions arc the conditions for an elec- 
tron system to be Galilean invariant. 
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The intcrband contribution can be simplified by noting that 

T bVl + ^ b - T a vl a = (nUk+ q ,6 - ^"k+qj + Oa< +q , a - T a ^J. (32) 

The second term on the right hand side of Eq. ( I32p can then be manipulated exactly as in the intra-band case to 
arrive at 

/W=^- Y\ / |(fo,k|^+Ao|a,k)| 2 A a (e F ,k)A b (e F ,k)^^ + ^nter (33) 

where 

*/W = ^- E / |(o,k- q|^ + A |6, k)| 2 A a (e F ,k - q)A t (6F,k) (T6Vk - b ~ TaVk ' a) ■ E . (34) 



When Galilean invariance is preserved the quasiparticle 
velocity and scattering times are the same for all bands, 
which implies that (5/3 = and hence that /Sinter = center- 
Although realistic materials arc not Galilean invariant, 
5/3 is nevertheless probably not significant because the 
term between parenthesis in Eq. (|34[) has an oscillatory 
behavior prone to cancellation. The degree of such can- 
cellation must ultimately be determined by realistic cal- 
culations for particular materials. 
With this proviso, we estimate that 







1 



8tts(] 



|(&,k|5 + A |a,k)| 2 ^ o (e F ,k)A 6 (e F ,k) 



q • Vd,a 

q v s 



(35) 



As long as 5(3 ~ is justified, the simplicity of Eq. 
in comparison to Eq. ([15]) or (JT7J) makes of the former the 
preferred starting point for electronic structure calcula- 
tions. Even when 5(3 ^ Q Eq. (|35[) may be an adequate 
platform for ab-initio studies on weakly disordered transi- 
tion metal fcrromagncts and strongly spin-orbit coupled 
ferromagnetic semiconductors* 2 ^ where (3 is largely de- 
termined by the intra-band contribution. Furthermore, 
a direct comparison between Eq. (|28|) and Eq. |35|) leads 
to the following observations. First, for nearly parabolic 
bands with nearly identical curvature, where the "drift 
velocity" is weakly dependent on momentum or the band 
index, we obtain (3 ~ (v^/v s )a and thus f3/a is roughly 
proportional to the ratio of the total spin density to the 
itinerant spin density, in concordance with predictions 
from toy models^ Second, if a/ (3 > for a system with 
purely electron- like carriers, then a/ [3 > for the same 
system with purely hole-like carriers because for a fixed 
carrier polarization v°[ and v s reverse their signs under 
m — > — m. However, if both hole-like and electron-like 
carriers coexist at the Fermi energy, then the integrand 
in Eq. ([35]) is positive for some values of a and negative 
for others. In such situation it is conceivable that a/ (3 be 
cither positive or negative. A negative value of (3 implies 
a decrease in magnetization damping due to an applied 
current. 
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FIG. 4: Comparison of a and /3 in (Ga,Mn)As for x — 0.08 
and p = 0Anm~ 3 . It follows that f3/a ~ 8, with a weak 
dependence on the scattering rate off impurities. If we use the 
torque correlation formula (Section VII), we obtain f)/a ~ 10. 



As an illustration of the foregoing discussion, in 
Fig. (gD we evaluate a/ (3 for (Ga,Mn)As. We find (3 to 
be about an order of magnitude larger than a, which is 
reasonable because (i) the local moment magnetization is 
larger than the valence band hole magnetization, and (ii) 
the spin-orbit coupling in the valence band decreases the 
transport spin polarization. Accordingly (3 is of the order 
of unity, in qualitative agreement with recent theoretical 
work2£. 



VII. TORQUE-CORRELATION FORMULA FOR 
THE NON-ADIABATIC STT 



Thus far we have evaluated non-adiabatic STT us- 
ing the bare vertex (a, k|5 + |6, k+q). In this section, 
we shall analyze an alternative matrix element denoted 
(a, k\K \b, k+q) (see below for an explicit expression), 
which may be better suited to realistic electronic struc- 
ture calculations. 16 ' 31 We begin by making the ap- 
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FIG. 5: M2DEG: comparing S and K matrix element ex- 
pressions for the non-adiabatic STT formula in the weakly 
spin-orbit coupled regime. Both formulations agree in the 
clean limit, where the intra-band contribution is dominant. 
In more disordered samples inter-band contributions become 
more visible and S and K begin to differ; the latter is known 
to be more accurate in the weakly spin-orbit coupled regime. 



FIG. 6: M2DEG: In the strongly spin-orbit coupled limit 
the intra-band contribution reigns over the inter-band contri- 
bution and accordingly S and K matrix element expressions 
display a good (excellent in this figure) agreement. Neverthe- 
less, this agreement does not guarantee quantitative reliabil- 
ity, because for strong spin-orbit interactions impurity vertex 
corrections may play an important role. 



proximation that the exchange splitting can be writ- 
ten as a constant spin-dependent shift H cx = AqS z . 
Then, the mean-field quasiparticlc Hamiltonian = 
H ex can be written as the sum of a spin- 



H 



kin 



H, 



(A') 



independent part , the exchange term, and the spin- 
orbit coupling Hso ■ 
the identity: 



With this approximation, we have 



(a,k|S+|6,k + q 
1 



A 
1 



[a,k\ 



[a,k\ 



|6,k + q) 
|6,k + q) 



(36) 



The last term in the right hand side of Eq. (j3"B"|) is the 
generalization of the torque matrix element used in ab- 
initio calculations of the Gilbert damping: 



(a,k|X|6,k + q) 



1 



(a,k| 



(37) 



Eq. (SB) implies that at q = (b, k\S + \a, k) ~ 
(b, k\K\a, k) provided that (E^ a — << A , which 

is trivially satisfied for intra-band transitions but less 
so for inter-band transitions^ For q^O the agreement 
between intra-band matrix elements is no longer obvi- 
ous and is affected by the momentum dependence of 
the band eigenstates. At any rate, Eq. (|29|) demon- 
strates that only q = matrix elements contribute to 
Antra! therefore Antra has the same value for S and K 
matrix elements. The disparity between the two formu- 
lations is restricted to Anter, and may be significant if 
the most prominent inter-band matrix elements connect 
states that are not close in energy. When they disagree, 




0.20 
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FIG. 7: GaMnAs: comparison between S and K matrix 
element expressions for Eli The disagreement between 
both formulations stems from inter-band transitions, which 
are less important as r increases. Little changes when E || z. 



it is generally unclear— whether S or K matrix elements 
will yield a better estimate of Antcr- The weak spin-orbit 
limit is a possible exception, in which the use of K ap- 
pears to offer a practical advantage over S. In this regime 
S generates a spurious inter-band contribution in the ab- 
sence of magnetic impurities (recall Section III) and it is 
only after the inclusion of the leading order vertex correc- 
tion that such deficiency gets remedied. In contrast, K 
vanishes identically in absence of spin-orbit interactions, 
thus bypassing the pertinent problem without having to 
introduce vertex corrections. 

Figs. ©- (0 display a quantitative comparison be- 
tween the non-adiabatic STT obtained from K and S, 
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both for the M2DEG and (Ga,Mn)As. Fig. (0) reflects 
the aforementioned overestimation of S in the inter-band 
dominated regime of weakly spin-orbit coupled fcrromag- 
nets. In the strong spin-orbit limit, where intra-band 
contributions dominate in the disorder range of interest, 
K and S agree fairly well (Figs. (JB]) and ((7])). Summing 
up, insofar as impurity vertex corrections play a minor 
role and the dominant contribution to (3 stems from intra- 
band transitions the torque-correlation formula will pro- 
vide a reliable estimate of f3. 



VIII. CONNECTION TO THE EFFECTIVE 
FIELD MODEL 

As explained in Section II we view the non-adiabatic 
STT as the change in magnetization damping due to a 
transport current. The present section is designed to 
complement that understanding from a different perspec- 
tive based on an effective field formulation, which pro- 
vides a simple physical interpretation for both intra-band 
and inter-band contributions to (3. 

An effective field TC cS may be expressed as the varia- 
tion of the system energy with respect to the magnetiza- 
tion direction Hf = -(l/s )dE/dn i . Here we approxi- 
mate the energy with the Kohn-Sham eigenvalue sum 



E 



^ ^k,a£k,a ■ 
k,a 



(38) 



The variation of this energy with respect to the magne- 
tization direction yields 



HI 



k,a 



de 



k . a 



an. 



dn Ka 

da,. 



"£k,a 



(39) 



It has previously been shown that, in the absence of cur- 
rent, the first term in the sum leads to intra-band Gilbert 
dampin g 15 ' 35 while the second term produces inter-band 
damping^! In the following, we generalize these results 



by allowing the flow of an electrical current, a and [3 may 
be extracted by identifying the the dissipativc part of the 
effective field with —adCt/dt — /3v s ■ Vf2 that appears in 
the LLS equation. 

Intra-band terms: We begin by recognizing that as the 
direction of magnetization changes in time, so does the 
shape of the Fermi surface, provided that there is an in- 
trinsic spin-orbit interaction. Consequently, empty (full) 
states appear below (above) the Fermi energy, giving rise 
to an out-of-cquilibrium quasiparticlc distribution. This 
configuration tends to relax back to equilibrium, but re- 
population requires a time r. Due to the time delay, 
the quasiparticle distribution lags behind the dynamical 
configuration of the Fermi surface, effectively creating a 
friction (damping) force on the magnetization. From a 
quantitative standpoint, the preceding discussion means 
that the quasiparticle energies e ka follow the magnetiza- 
tion adiabatically, whereas the occupation numbers n ki£l 
deviate from the instantaneous equilibrium distribution 
h, a via 



n k,a — fk,a — 7k : a 



dfi 



dfu 



dt 



Or 



dfk,g 

<9k 



(40) 

where we have used the relaxation time approximation. 
As we explain below, the last two terms in Eq. ([40]) do 
not contribute to damping in the absence of an electric 
field and have thus been ignored by prior applications of 
the breathing Fermi surface model, which concentrate on 
Gilbert damping. It is customary to associate intra-band 
magnetization damping with the torque exerted by the 
part of the effective field 



H 



eff 



s ° k 

k,a 



de^a 
80. 



(41) 



that is lagging behind the instantaneous magnetization. 
Plugging Eq. (gDJ) in Eq. (gl|) we obtain 



J 



eff 

intra, i 



so 



E 

k . a 



dtk,a , df k . a <9e k . a <9e k . a dflj 
dfli <9ek,a dQ,i dilj dt 



+ T a f 



dfk,g dt^g dt^g 8Qj 



5/k,a dek,a de k , a 



r 



(42) 



where a sum is implied over repeated Latin indices. The 
first term in Eq. ()42|) is a contribution to the anisotropy 
field; it evolves in synchrony with the dynamical Fermi 
surface and is thus the reactive component of the effective 
field. The remaining terms, which describe the time lag 
of the effective field due to a nonzero relaxation time, are 
responsible for intra-band damping. The last term van- 
ishes in crystals with inversion symmetry because k = eE 
and 9e/<9k is an odd function of momentum. Similarly, 



if we take r = <9e(k)/9k the second to last term ought to 
vanish as well. This leaves us with the first two terms in 
Eq. ( I42p , which capture the intra-band Gilbert damping 
but not the non-adiabatic STT. This is not surprising as 
the latter involves the coupled response to spatial varia- 
tions of magnetization and a weak electric field, render- 
ing linear order in perturbation theory insufficient (see 
Appendix A). In order to account for the relevant non- 
linearity we use r = <9e(k — ev-Er)/9k in Eq.( l4"2"|) . where 
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de(k) 1 <9k. The dissipative part of Tif^ tra then reads 



We thus concentrate on the second term in Eq. ([39)) 



H 



eff,damp 

intra, i 



— V T k , a 
k,a 



9ek a dVli d£lj 



<9fL 



dt A ' a dr l 
(43) 



1 dn Ka 

«n ^ ^ 



k,a 



(46) 



where v i da = eTa^i^ 1 ) 1 ^^ is the "drift velocity" cor- Multiplying Eq. (|46|) with dti/dt we get 
responding to band a. Eq. (|43[) may now be identified 
with — ai n tra<9r2/dt — /?intraVs • Vf2 that appears in the 
LLS equation. For an isotropic system this results in 



-j_/cff ,damp 
Writer ' u t 



■ an 



Centra — — 



A, 



so 
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so 



Tk.a 



> Tk,o -q 
k.a.z 



d/k,a / de k , Q 
<9e k , a V dtli 

dfk, a ( de k , a 

on, 



(44) 



Since {{S x ,H so }) = a^(exp(iS S! 0) J ff flO exp(-iS a: ^)) = 
de/d(j) for an infinitesimal angle of rotation <^> around 
the instantaneous magnetization, /3 in Eq. ([4"4"|) may be 
rewritten as 

Antra = ^ £ T ^ ^ I (k, a|tf|k, a) | 2 (45) 



k.a 



9e k .a 



q • v s 



where K = [S + , H so ]/Aq is the spin-torque operator in- 
troduced in Eq. ( I37[) and we have claimed spin rota- 
tional invariance via | ([S x , H so \) \ 2 = \([S y , H so })\ 2 . Using 
df/de ~ — 8{e — e^-) and recalling from Section VII that 
K a ,a = S+ a , Eq. IJ25J) is equivalent to Eq. ([301); notc tnat 
the product of spectral functions in the latter yields a 
factor of 4-7TT upon momentum integration. These obser- 
vations prove that Antra describes the contribution from 
a transport current to the "breathing Fermi surface" type 
of damping. Furthermore, Eq. (j44| highlights the impor- 
tance of the ratio between the two characteristic veloci- 
ties of a current carrying fcrromagnct, namely v s and v c i- 
As explained in Section VI these two velocities coincide 
in models with Galilean invariance. Only in these arti- 
ficial models, which never apply to real materials, does 
a = {) hold. 

Inter-band terms: The Kohn-Sham orbitals are effec- 
tive eigenstates of a mean-field Hamiltonian where the 
spins arc aligned in the equilibrium direction. As spins 
precess in response to external rf fields and dc trans- 
port currents, the time-dependent part of the mean-field 
Hamiltonian drives transitions between the ground-state 
Kohn-Sham orbitals. These processes lead to the second 
term in the effective field and produce the inter-band 
contribution to damping. 



so 



£k,a dn a ^/dh ■ dCl/dt 



= -— y^ek.adn aM /dt. 

k,a 



(47) 



The rate of change of the populations of the Kohn- 
Sham states can be approximated by the following master 
equation 



dt 



= Wa,b(«k,a - «-k',b), 



(48) 



h,k' 



where 



W a , b = 2TT\(b 1 k , \A S x \a,k)\ 2 5 k ' M+cl S(e bM/ - e aM - uj) 

(49) 

is the spin-flip inter-band transition probability as dic- 
tated by Fermi's golden rule. Eqs. (gSJ) and l[4"9"]) 
rely on the principle of microscopic reversibility^ and 
are rather ad hoc because they circumvent a rigorous 
analysis of the quasiparticle-magnon scattering, which 
would for instance require keeping track of magnon occu- 
pation number. Furthermore, quasiparticle-phonon and 
quasiparticle-impurity scattering are allowed for simply 
by broadening the Kohn-Sham eigenenergies (see below). 
The right hand side of Eq. (|4"B"]) is now closely related 
to inter-band magnetization damping because it agrees^ 
with the net decay rate of magnons into particle-hole 
excitations, where the particle and hole are in different 
bands. Combining Eq. (|4"7| and (|4"8"]l and rearranging 
terms we arrive at 



Winter ' ~ 



u k,k',a.b 



W a ,b(nk,a - nk',&)(£k,a - e k',o)- 

(50) 

For the derivation of center it is sufficient to approximate 
Ti\^ a as a Fermi distribution in Eq. (|50[k here we ac- 
count for a transport current by shifting the Fermi seas as 
«k,a -* "k.a - ev k , a ■ ET k: a<9n k ,a/<9ek,a, which to leading 
order yields 
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|(6,k + q|A 5 + |a, k)| 2 (5(e b ,k+ q -ea,k-w) 



k,a,6 



<9"-k,a 
<9ek,a 



(— w + eT4 j0 ) 



Sttsq 



J2 I (6, k + q| Ao^ + |a, k) | 2 A Q (k, e F )A(k + q, e*.)(-w + eH.a) 



(51) 



k,a,6 



where we have used S x = (S + + S~)/2 and defined VJ,. a = evk +qj £, • Erk+ q ,6 — evj^a • Etj. a . In the second line of 
Ea.f I5ip we have assumed low temperatures, and have introduced a finite quasiparticle lifetime by broadening the 
spectral functions of the Bloch states into Lorentzians with the convention outlined in Eq. (|14[) : 6(x) — ► A(x)/(2tv). 
Identifying Eq.([5l]) with (—ai nteT d t Q - /3 in tcr(v s • V)fi) • d t Cl = -a- lntcr Lu 2 + /3- mteT cj(q ■ v s ) we arrive at 



Winter 



0n 



— 12 \(b,k + q\A S + \a,k)\ 2 A a (k,e F )A b (k + q,e F ) 



8tTSq 



l 



87rs q • v s 



]T J2 \ (b,k + q\A S+\a,k)\ 2 A a (k,e F )A b (k + q,e F )V b , a 



(52) 



a.b^a k,a,b 

in agreement with our results of Section II. 



r 



fers a physically sensible result in the weak spin-orbit 
limit but its quantitative reliability is questionable un- 
less the prominent inter-band transitions connect states 
that are close in energy. Finally, we have extended the 
breathing Fermi surface model for the Gilbert damping 
to current carrying ferromagnets and have accordingly 
found a complementary physical interpretation for the 
intra-band contribution to /?; similarly, we have applied 
the master equation in order to offer an alternative inter- 
pretation for the inter-band contribution to /3. Possible 
avenues for future research consist of carefully analyzing 
the importance of higher order vertex corrections in /?, 
better understanding the disparities between the differ- 
ent approaches to v s , and finding real materials where 
a/P is negative. 



IX. SUMMARY AND CONCLUSIONS 

Starting from the Gilbert damping a and including the 
influence of an electric field in the transport orbitals semi- 
classically, we have proposed a concise formula for the 
non-adiabatic spin transfer torque coefficient [3 that can 
be applied to real materials with arbitrary band struc- 
tures. Our formula for (3 reproduces results obtained 
by more rigorous non-linear response theory calculations 
when applied to simple toy models. By applying this ex- 
pression to a two-dimensional electron-gas ferromagnct 
with Rashba spin-orbit interaction, we have found that 
it implies a conductivity-like contribution to /3, related to 
the corresponding contribution to the Gilbert damping a, 
which is proportional to scattering time rather than scat- 
tering rate and arises from intra-band transitions. Our 
subsequent calculations using a four-band model have 
shown that intra-band contributions dominate in ferro- 
magnetic semiconductors such as (Ga,Mn)As. We have 
then discussed the a/f3 ratio in realistic materials and 
have confirmed trends expected from toy models, in ad- 
dition to suggesting that a and P can have the oppo- 
site sign in systems where both hole-like and electron-like 
bands coexist at the Fermi surface. Afterwards, we have 
analyzed the spin-torque formalism suitable to ab-initio 
calculations, and have concluded that it may provide a 
reliable estimate of the intra-band contribution to /?; for 
the inter-band contribution the spin-torque formula of- 

I 

APPENDIX A: QUADRATIC SPIN RESPONSE TO AN ELECTRIC AND MAGNETIC FIELD 

Consider a system that is perturbed from equilibrium by a time-dependent perturbation V(t). The change in the 
expectation value of an operator 0(t) under the influence of V(t) can be formally expressed as 
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S(0(t)) - (* \rf(t)O(t)U(t)\* ) - (*o|0(f)|tt ) 



(Al) 
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where |\&o) is the unperturbed state of the system, 

U(t) = Texp 



V(t')dt' 



(A2) 



is the time-evolution operator in the interaction representation and T stands for time ordering. Expanding the 
exponentials up to second order in V we arrive at 



dt>([0(t),V(t')})-± I 



For the present work, 0(t) — > S a (a — x, y, z) and 



V(t) 



J dry A(r,t) + J drS ■ H ext (r,t), 



(A3) 



(A4) 



where A is the vector potential, 7i cx t is the external magnetic field, and j is the current operator. Plugging Eq. (|A4I 
into Eq. (|A3[) and neglecting 0(A 2 ), 0(H^ xt ) terms we obtain 



SS°(x) =J2 J dx'x a si Ab (*')+Y, J dx' x a s b s nU^) + Y, j dx'dx" X a ^ c .A\x')Hl 



(A5) 



where x = (r, t) and J dx' = dt' J dv' . The linear and quadratic response functions introduced above are defined 



as 



xt 3 M = i{[s a (x)j b (x')]e(t-t') 

Xs b s (x,x') = i([S a (x),S b (x')]Q(t-t') 
Xfs^x'^x") = ([[S a (x)J b (x')] ,S c (x")] Q(t-t')e(f-t") 

+([[s a {x),s b {x")] ,f(x')] ®(t - t")e(t" - 1) 



(A6) 



where we have used T [F(t)G(t')} = F{t')G(t")®{t' -t") + G(t")F(t')Q{t" -t'), 9 being the step function. X s,j is the 
spin density induced by an electric field in a uniform ferromagnet, and it vanishes unless there is intrinsic spin-orbit 
interaction, xs.s is the spin density induced by an external magnetic field. Xs,S,j is the spin density induced by the 
combined action of an electric and magnetic field (sec Fig. {8} for a diagrammatic representation); this quantity is 
closely related to (v s • q)x^ 2 \ introduced in Section II. 



APPENDIX B: FIRST ORDER IMPURITY VERTEX CORRECTION 



The aim of this Appendix is to describe the derivation of Eq. ([?!"]) . We shall begin by evaluating the leading order 
vertex correction to the Gilbert damping. From there, we shall obtain the counterpart quantity for the non-adiabatic 
STT by shifting the Fermi occupation factors to first order in the electric field. The analytical expression for the 
transverse spin response with one vertex correction is (see Fig. 

xf- {1) = - V ^ T Y] I u l G a (iu; n ,k)S+ b G b (iu; n + iLJ : k+q)Sl b ,G bl (iLJ n + iu J ,k' + q)S^ a ,G a ,(itJ n X)Sl a . (Bl) 

C3 

r v.a 



FIG. 8: Feynman diagram for \s,s,j- The dashed lines correspond to magnons, whereas the wavy line represents a photon. 
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FIG. 9: Feynman diagram for the first order vertex correction. The dotted line with a cross represents the particle-hole 
correlation mediated by impurity scattering. 



where V is the volume of the system and the minus sign originates from fermionic statistics. Using the Lchmannn 
representation of the Green's functions G and performing the Matsubara sum we get 



~QP,(i) 
X+ - 



A 2 



u l 2 Rc 



k,k' 



Q+ QX Q— Ql 

J a,b°b,b' J &'.a"~V,a 



xA b (e 2 ,k + q)A bl (e' 2 ,k' + q) 



/(ei) 



(ei - e[){iLu + ei - e 2 ){iu) + e x - e 2 ) 



■2- 



(B2) 



where twice the real part arose after absorbing two of the terms coming from the Matsubara sum. Next, we apply 
iuj — > lu + i0 + and take the imaginary part: 



~QP,U) = 



u l Re 



k,k' 



de^deide, ^ ^ (4, k ' )A (e 2 , k + q)A 6 , (e 2 , k' + q) 

(B3) 



, /(ei) 



<5(w + ei — £2) 5(u + t\ — e' 2 ) 



lo + ei - e 2 



u + t\ - e 2 



(2tt)< 

LU — > — W, 

q-> -q 



where we used l/(x — ir/) = PV(l/x) + inS(x), and invoked spin-rotational invariance to claim that terms with 
S* b Sl b ,Sf, a ,S l a , will vanish. Integrating the delta functions we arrive at 



A 2 r 
z Jk,k' 



de'^de'z f(e 2 )A a (e 2 ,k)A a , (ei, k') 



(2tt) 3 



(e 2 -e 2 )( e2 -6' 1 ) 



[A h (e 2 + w, k + q)A b , (e 2 + cj, k' + q) + A&(e' 2 + w, k + q)A b > (e 2 + u>, k' + q)] - 



(B4) 



The next step is to do the and e' 2 integrals, taking advantage of the fact that for weak disorder the spectral 
functions are sharply peaked Lorcntzians ( in fact at the present order of approximation one can take regard them as 
Dirac delta functions). The result reads 



~QP,(i) 
X+,- 



a 2 r 

/ tfRe I 

1 Jk,k' 



de 2 f(e 2 )A a {e 2l k) 
-00 27r e 2 - e k ' „/ 



A b (e 2 +u,k + q) A 6 ,(e 2 +u;,k' + q) 



e 2 + LU — £k'+q^' 



£ 2 + LU — £k+q,6 



- (w -> -u, q -> -q) 
By making further changes of variables, this equation can be rewritten as 

~A b (e 2 +w,k + q) A b ,{e 2 + lu, k' + q) 



~QP,(1) = F A 



tt <Re[...] r &2 (/(e2)-/(e2+^))^ a (e2,k) 
2 ./k.k' J-00 27T e 2 - e k ' a / 



£ 2 + W — Ek'+q,b' 



£ 2 + W — £k+q,6 



(B5) 



(B6) 



This is the first order vertex correction for the Gilbert damping. In order to obtain an analogous correction for the 
non-adiabatic STT, it suffices to shift the Fermi factors in Eq. (|B6[) as indicated in the main text. This immediately 
results in Eq. (|2Tj) . 



APPENDIX C: DERIVATION OF EQ. (^6)) 

Let us first focus on the first term of Eq. (fTT)) . namely 

E iQj J [I (a, k|5+|6, k) | 2 + I (a, k|5" \b, k) | 2 ] A a 4,<X >6 r k , 



(CI) 
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We shall start with the azimuthal integral. It is easy to show that the entire angle dependence comes from v % tfi oc kikj, 
from which the azimuthal integral vanishes unless i = j. 

Regarding the \k\ integral, we assume that A/cf,Ao,1/t << ep; otherwise the analytical calculation is complicated 
and must be tackled numerically. Such assumption allows us to use J k — > N2D J_ OD de. For inter-band transitions 
(a 7^ b), A a A' b contributes mainly thru the pole at £f,o, thus all the slowly varying factors in the integrand may be 
set at the Fermi energy. For intra-band transitions (a = b), A a A' a has no peak at the Fermi energy; hence it is best 
to keep the slowly varying factors inside the integrand. 

The above observations lead straightforwardly to the following result: 



EiQi J [| (a, k\S+\b, k) | 2 + I (a, k\S-\b, k) | 2 ] A a A' b vl a vi b r Ka 

( e F,-T_r + - e F , + T + r_) 



E 



E 



8m_i_m_ 



A 2 
b 2 



b 3 



i 2 1 A 2 fc 2 



2 b 2 



A 2 
b 2 



lX 2 k 



2 1.2 



2 b 2 



A 2 
b 2 



(C2) 



The second and third line in Eq. (|C2[) come from inter-band and intra-band transitions, respectively. The latter 
vanishes in absence of spin-orbit interaction, leading to a 2D version of Eq. (j20|) . Since the band-splitting is much 
smaller than the Fermi energy, one can further simplify the above equation via r+ ~ r_ — > r. 
Let us now move on the second term of Eq. (|17[) . namely 



EiqjJ-Re [(b,k\S'\a,k)(a,k\S + d kj \b,k) + (S+ 4-> S~j\ A a A b v l Ka T^ a 



(C3) 



Most of the observations made above apply for this case as well. For instance, the azimuthal integral vanishes 
unless i = j. This follows from a careful evaluation of the derivatives of the eigenstates with respect to momentum; 
dk = sin(6>) cos(9)kj /k 2 (0 < 9 < tt/2) is a useful relation in this regards, while d k 4> plays no role. As for the |fc| 
integral, we no longer have the derivative of a spectral function, but rather a product of two spectral functions; the 
resulting integrals may be easily evaluated using the method of residues. The final result reads 

E t q [ Re[(b,k\S-\a,k)(a,k\S + d k] \b,k) + (S+ ^ S-)] A a A b vl a r Ka 
Jk 

~ -E- q 
+ Eq 



ni 



32m_ 
m 



A 2 fcp A 2 



1 + — 



Am j 



\2i,2 a 2 
A ftp^o 2 
7~i 



If 1 



m 
32m. 
m A 2 fc|A^ 2 



A 2 fc 2 A 2 
b 6 



l + T -± 



4m_ 



fr 4 



(C4) 



The first line in Eq. (|C4[) stems from inter-band transitions, whereas the second comes from intra-band transitions; 
both vanish in absence of SO. Once again we can take r + ~ r_ — > r. Combining Eqs. (|C2[) and (|C4p one can 
immediately reach Eq. (|26p . 
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